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Modification of surface and bulk properties of solids by irradiation with ion beams is a widely used 
technique with many applications in material science. In this study, we show that nano-hillocks on C3F2 
crystal surfaces can be formed by individual impact of medium energy (3 and 5 MeV) highly charged ions 
(Xe^^^ to Xe^^^) as well as swift (kinetic energies between 12 and 58 MeV) heavy xenon ions. For very slow 
highly charged ions the appearance of hillocks is known to be linked to a threshold in potential energy (Ep) 
while for swift heavy ions a minimum electronic energy loss per unit length (Sg) is necessary. With our 
results we bridge the gap between these two extreme cases and demonstrate, that with increasing energy 
deposition via Sg the Ep-threshold for hillock production can be lowered substantially. Surprisingly, both 
mechanisms of energy deposition in the target surface seem to contribute in an additive way, which can be 
visualized in a phase diagram. We show that the inelastic thermal spike model, originally developed to 
describe such material modifications for swift heavy ions, can be extended to the case where both kinetic and 
potential energies are deposited into the surface. 



The interaction of energetic ions with solid surfaces has been a center of attention in material science over the 
last century, not least because profound knowledge of the underlying physical mechanisms has contributed 
greatly to technological applications. A number of analytical methods, e.g. Rutherford backscattering 
spectrometry (RBS) or secondary ion mass spectrometry (SIMS), could not have been developed without this 
continuous research. Energy deposition by energetic ions in matter also has important applications in other fields 
like astrobiology (synthesis of amino acids), nuclear fusion (wall erosion in fusion experiments), space explora- 
tion (radiation hardness of electronic devices) and others. 

Great interest has been drawn by the use of ions as possible tools for lithographic modifications or "nanos- 
tructuring" of surfaces. Due to technical restrictions, earlier studies in the field of nanostructuring were conducted 
with ions at relatively high kinetic energies and low ion charge states, where interaction processes are governed 
primarily by the ion mass and kinetic energy. Investigations have shown that swift heavy ions (SHI) with kinetic 
energies in the MeV-GeV range can induce severe structural modifications at the surface^ and in the bulk^. These 
modifications not only include the formation of surface nanostructures, such as protrusions^ or craters^, but also 
phase transitions by swift heavy ions (e.g., from crystalline to amorphous^ or crystalline to crystalline^) and the 
creation of latent tracks, which can reach several tens of micrometers deep into the solid^. 

The interaction of individual slow highly charged ions (HCI) in the keV energy regime with solid surfaces can 
also induce surface modifications on a nanometric scale^ which appear as either hillocks, pits or craters^"^^. These 
modifications result from the deposition of potential energy carried by slow HCI (i.e. the sum of the binding 
energies of all missing electrons) into the target electronic system. Within a few femtoseconds (fs), electrons from 
the target are captured into highly excited states of the projectile^^. The energy (several keV/nm^) deposited into 
the electron target by emission of fast electrons is not precisely known^^. Similarly, hillocks created by SHI with 
impact energies of some hundreds of MeV appear at the surface for an energy density of typically 10 keV/nm^. 
These energies provided by ion-electron collisions in the fs time range^^"^^, have been measured^^"^^. In this case, 
fast electrons are created by direct collisions with the ions and their energy distribution can be approximated by an 
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analytical formula proposed by Waligorski et al.^^ that fits the results 
of Monte Carlo calculations^^. We conclude that independent of the 
energy deposition channel to the electrons, i.e. potential energy (Ep) 
deposition by HCI or electronic energy loss per length unit (Sg) by 
SHI, the target electronic system is highly excited in a first step by the 
generation of fast electrons. The electron energy relaxation towards 
the atoms occurs in a second step, which suggests that the material 
modifications induced at the surface have a common underlying 
mechanism, as pointed out by Aumayr et al.^. Interpreting their 
results for CaF2, El-Said et al.^ were the first to put forward that 
hillock formation by slow HCI may be described by the inelastic 
thermal spike modeP^ originally developed to predict material mod- 
ifications induced by SHI. Along the lines of this model, the sharp 
threshold of potential energy for hillock formation by slow HCI has 
been linked to a solid-liquid phase transition (nano - melting) 

In this paper we firstly establish the so far missing quantitative link 
between HCI and SHI induced nano-hillocks. By bombarding CaF2 
single crystal surfaces with medium energy (3 and 5 MeV) multiply 
charged xenon ions which deposit a considerable amount of both 
potential and kinetic energy into the surface, we can demonstrate 
that both mechanisms of energy deposition into the target surface are 
found to contribute in an additive way. Based on this relation, the 
inelastic thermal spike model, originally developed to describe such 
material modifications for swift heavy ions, can be extended to 
nanostructure formation by highly charged ions. 

Results 

Hillock formation and surface roughening by medium energy 
HCI. Fig. 1 shows AFM images of samples irradiated by 5 MeV 
Xe'^+ (a), Xe''+ (b) and Xe''+ (c) ions. For Xe'^+ (Fig. la), no 
hillocks are observable and the irradiated sample stays flat with a 
mean roughness of 0.07 nm. The same is true for surfaces irradiated 
by ions in lower charge states (Xe^°^, Xe^^^). However, nano-hillocks 
clearly appear on surfaces after irradiation by Xe^^^ (Fig. lb), Xe^^^ 
and Xe^°^. For these cases, the number of hillocks per area corres- 



ponds well to the fluence of incident ions within the experimental 
error. Samples irradiated by 3 MeV Xe^^^ ions (and ions in lower 
charge states) do not exhibit hillock formation, while for Xe^^^, 
hillocks appear clearly but with the number of hillocks reduced to 
1.5% of the incoming ion fluence. When using 3 MeV Xe^°^ ions as 
projectiles, hillocks are clearly visible and their number is equal to the 
ion fluence. 

Fig. Id shows the roughness of the irradiated samples, as calcu- 
lated by the Gwyddion code^^, plotted against the charge state of the 
Xe ions for different kinetic energies (3 and 5 MeV, plus one irra- 
diation performed at 0.54 MeV). This illustrates a surprisingly sharp 
but impact energy dependent transition from flat to nanostructured 
CaF2 surfaces. The red dotted line which is just a link between Xe^^^ 
and Xe^°^ for 3 MeV irradiation is to guide the eye. Right at the 
threshold (Xe^^+ for 5 MeV and Xe^^+ for 3 MeV), the surface rough- 
ness increases significantly to values 4-10 times larger than the non- 
irradiated sample. Summing up the results it can be said that for Xe 
ions of 5 MeV, the threshold for nanostructuring lies between charge 
states 21 + and 22 + , which corresponds to a potential energy of Ep5 
— 5.5 keV^^. For lower kinetic energies of 3 MeV, the threshold is 
given by Xe^^^ and Xe^^^ or a potential energy of Ep3 ~ 8.5 keV^^. 
Both threshold values of Ep in this energy regime are considerably 
smaller than the ones in the keV regime as reported by El- Said et al.^ 
(between Xe^^^ and Xe^^^, corresponding to Epo ~ 11.2 keV). 

With increasing kinetic energy, nuclear and electronic energy loss 
per unit length of the projectile ion evolve differently. For energies 
between 0.004^ and 5 MeV, the nuclear energy loss is nearly constant, 
varying only from 1.15 keV/nm to 1 keV/nm respectively^^. Such a 
15% decrease of the nuclear energy loss within increasing kinetic 
energy cannot account for the observed decrease of the Ep threshold 
by 50% from 11.2^ to 5.5 keV. But in the same energy range, Sg 
increases considerably^°'^\ Since the electronic system of the target 
is excited by Sg, independent of whether the energy is supplied by SHI 
or HCI, we conclude that a decrease of the Ep threshold has to be 
counter-balanced by an increased energy deposition via Sg. 
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Figure 1 | CaF2 surfaces irradiated by 5 MeV Xe ions with a fluence of 500 ions/^im^. The AFM images show an area of 1 X 1 [im^. A surface irradiated 
by (a) Xe^^^with no hillocks visible, (b) Xe^^^ with —520 hillocks, (c) Xe^°^ with —410 hillocks, (d) presents the evolution of the surface roughness 
versus charge state for samples before and after irradiation with 5, 3 and 0.54 MeV Xe ions. The lines are there to guide the eye. 
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Hillock formation by SHI and the velocity effect. To determine the 
Se threshold for hillock creation by SHI, targets of CaF2 were 
irradiated with some tens of MeV of Xe ions delivered by the 
IRRSUD beam line of the GANIL facility^l Al foil degraders were 
used to reduce the kinetic energy of the initial 92 MeV Xe ion beam 
to values ranging between 12 and 58 MeV. The fluences applied were 
5 X 10^ Xe/cm' and 1.5 X W Xe/cml Using AFM in tapping 
mode^^, the surface of each sample was analyzed after irradiation 
to extract the number of hillocks per cm^. The efficiency for hillock 
production (number density of hillocks divided by incident ion 
fluence) was deduced for irradiations at two fluences. Their mean 
values versus Sg are plotted in Fig. 2. The Sg values were calculated 
using the CasP code^^ since its predictions of Se values are more 
realistic as compared to the ones by SRIM^^ in the low MeV energy 
regime. 

By extrapolation to zero efficiency, a linear fit, as previously done 
like in the GeV energy regime^^, was applied to the three points at Sg 
= 3.6, 5.9 and 9.5 keV/nm, leading to Se threshold of Set =2.75 ± 
0.40 keV/nm. This threshold value is by a factor of ~2 smaller than 
the Se threshold (~5 keV/nm) previously reported for SHI in the 
GeV regime^^. Such a difference in the Se threshold can in principle 
result from the so-called "velocity effect"^°: For SHI impact on a 
surface with different kinetic energies but with the same Se value, 
the energy deposited into the electronic system at lower ion velocity 
is transferred to the lattice atoms more efficiently than at higher 
energies. From our results we infer that CaF2 is sensitive to the 
velocity effect^ ^ (as questioned by Szenes^^). 

In Fig. 3 we have summarized the results of our investigations in a 
"phase diagram" with Se and Ep as state variables. Cases where hil- 
locks have been found after irradiation are displayed by full symbols 
(region B), while open symbols represent irradiations where the CaF2 
surface stayed flat and showed no effect of the irradiation (region A). 
Additional data from slow HCI impact^ have been included in this 
plot. The Se values (-0.01 keV/nm for 3 keV, 0.72 keV/nm for 
3 MeV and 1.46 keV/nm for 5 MeV Xe ions) were calculated with 
the CasP code^\ The SHI threshold Set of 2.75 keV/nm (semi-black 
square in Fig. 3) for hillock formation is also given in the phase 
diagram, assuming a potential threshold of 0.62 keV corresponding 
to the expected equilibrium charge state of the Xe at 9 MeV^\ The 
borderline between region A (no hillocks) and region B (hillocks) in 
the Se-Ep diagram (Fig. 3) is a linear fit of the Ep thresholds and the 
negative slope points out an additive contribution of Se and Ep to 
nano -hillock formation on the surface of CaF2. A decrease of Ep, 
transferred to the target clearly can be compensated by an increase 
of Se and vice versa. The extrapolation of this line to Ep = 0 gives Seo 
of 2.75 ± 0.40 keV/nm (semi-black square in Fig. 3) is in agreement 
with the Set as deduced from SHI experiments. 
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Figure 2 | Mean efficiency for hillock formation by SHI versus Sg from 
the present experiment with Xe projectile energies in the MeV energy 
regime (-0.2 MeV/u). 
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Figure 3 | Phase diagram of hillock formation on CaF2 irradiated by Xe 
ions with electronic energy loss and potential energy as state variables. 

The corresponding kinetic energies and Xe charge states are given on the 
opposite axes, respectively. Open symbols represent cases where no 
hillocks have been found (region A), while full symbols mean the 
appearance of hillocks (region B). The blue diamonds are previous results'^. 
The red circles stand for the present experiments. The semi-black square is 
the Sg threshold deduced from SHI irradiations. 

Data analysis: link between HCI and SHI at the thresholds of 
hillock formation. The additivity between Sg and Ep derived from 
Fig. 3 can be expressed in a quantitative formula taking into account 
that only a fraction F^^ of Ep is deposited in a depth d near the 
surface^^. 



Se + (F/d) X Ep = const. = Seo = Set = 2.75 keV/nm 



(1) 



F can either be derived from calorimetric measurements^^ or by 
measuring the energy carried away by emitted electrons^^ From 
our results (equ. 1), a ratio of the depth d in nanometer to the 
fraction F can be derived for CaF2 in beam energy range between 
0.004 and 9 MeV: 



d/F = 4.1 + 0.3 nm = Epo/Seo 



(2) 



where Epo is the potential threshold at Sg = 0 and Sgo (= Set) the 
electronic energy loss for Ep = 0. 

Inelastic thermal spike (i-TS) description. In the following we 
present an extension of the inelastic thermal spike (i-TS) modeP^, 
originally developed to describe material modifications induced by Se 
delivered by SHI, to medium and low energy HCI, where potential 
energy is deposited. As a central assumption, the threshold for hillock 
formation is linked to a solid-Uquid phase transition (nano- 
melting)^'^ inside the target material. 

i-TS model simulation for SHI. The original i-TS model describes a 
transient thermal process based on heat transport equations that 
govern the heat diffusion in time and space (radial distance from 
the ion path) in the electronic and atomic subsystem and their 
coupling by the electron-phonon constant g. For SHI the two 
equations are solved numerically in cylindrical geometry. The 
initial energy distribution of the electrons is derived from Monte 
Carlo calculations and implemented as an analytical formula^^. 
The electron-phonon coupling constant g is linked to the electron- 
phonon mean free path that defines the mean length of energy 
diffusion of the electrons before its transfer to the atoms. The model 
was successfully applied for SHI in the GeV energy regime for several 
solids, independent of whether they are metals or insulators^^ ''"'^ 
references therein^ dlrectty applied to CaF2'' wlth X equal to 
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Table 1 | Evolution of the radius of the molten phase at the surface 
(Rm) and the melting depth (dm) for different conditions of the initial 
radius (Rp) and depth (d) of the energy (F*Ep) transferred to the 
electrons 
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3.8 nm^'^ It is assumed that tracks and hillocks appear, if an 

energy of = 0.58 eV/atom, which is necessary to melt the 
material, is surpassed. Such calculated Sg threshold for material 
modification is predicted by the i-TS model to be 5 keV/nm in the 
GeV regime in close agreement with experiments^^. Applying the i- 
TS model to heavy ions in the MeV range^° (beam energy of 
0.08 MeV/u, leading to an initial energy transfer to the electrons of 
the target in a cylinder radius of 1 nm^^), a lower threshold of 2.6 ± 
0.4 keV/nm is calculated for hillock formation, which is in excellent 
agreement with our experimental results presented in Fig. 2. The 
"theoretical" error of 0.4 keV/nm basically reflects the 
uncertainties in the initial energy distribution of the electrons. 

Extension of the i-TS model for HCI - a 3D development. In order 
to take into account that the energy transferred to the electrons by Ep 
and Se is not constant along the ion path, the 3D i-TS model 
developed for metallic materials^^ has been adapted for insulators, 
following the assumptions made for the electron subsystem as 
described in the previous section^^ Moreover it is presumed that 
the electron-phonon mean free path (k = 3.8 nm) is the one 
previously used for the description of the track induced by SHI. It 
is also assumed that a fraction (F) of Ep is homogeneously deposited 
into a cylindrical volume below the surface, which is characterized by 
a depth (d) and a radius (Rp). These assumptions are used to simplify 
the first application of the numerical solution of the i-TS model to 
describe the threshold of hillock formation by HCI. Moreover, in a 
layer of depth d below the surface, Sg and F*Ep are deposited 
additively, while for larger depth Sg is the only source of energy. 
The radial distribution of the potential energy has been estimated 
by Lemell et al.^^, showing that the most efficiently interacting 
electrons are the ones with energies lower than 400 eV. These 
electrons are confined in a cylinder radius Rp between 0.5 and 
2 nm. We can see from this that in order to perform the 
calculations, it is necessary to estimate d, Rp and F. 

Several calculations were made to investigate the influence of Rp 
and d on the determination of the radius of the molten phase (R^) at 
the surface and on the melting depth (d^) for a Xe ion energy of 
5 MeV (Se =1.4 keV/nm), see Table 1 and Fig. 4. For the same d 
value (3 nm) and the same deposited potential energy (F*Ep = 
8.2 keV), an increase of Rp from 0.6 to 1.6 nm causes a 10% and 
20% decrease of R^ and d^, respectively (see Fig. 4a-c). Comparing 
Fig.4b and 4d for the same F*Ep (8.2 keV) and Rp (1 nm) values, a 
decrease of the depth from 3 to 1 nm induces a decrease of R^ by 4% 
and dm by 20%, respectively. The small variation of R^ and d^ results 
from the smearing out of the energy transferred to the electron sub- 
system of the target by the electron -electron interaction prior its 
transfer to the atoms due to the fact that Rp and d are smaller than 
the electron-phonon mean free path value. 

We will now use the i-TS model to describe the threshold of hillock 
appearance (Fig. 3) taking into account that the ratio d/F is constant 
in the range of potential energy and electronic energy loss investi- 
gated here. Fig. 5a presents the resulting maximum energy E^ trans- 
ferred to a target atom for 3 MeV Xe ion impact on CaF2 versus the 
deposited Ep in different depths. The choice of d values was done 
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Figure 4 | Effect of different values of Rp and d of the cylinder (dashed 
blue rectangle) where the energy is transferred to the electrons compared 
to the calculated volume of the molten matter (the red line characterizes 
the evolution of the melting depth versus the melting radius). The 

irradiation parameters are Xe ions with 5 MeV, leading to Sg = 1.4 keV/ 
nm and a deposited energy of F*Ep = 8.2 keV. 

according to the predictions of Lemell et al.^^ and to the measurement 
of non- equilibrium energy loss^^ in CaF2. The deposited potential 
energies needed to melt CaF2 are deduced when E^ reaches E^ 
(0.58 eV/at) and are equal to 3. 1,4.1, 5.3 and 6.9 keV for the different 
depths (1, 2, 3 and 4 nm respectively). The respective values of the d/ 
F ratio divided by the corresponding potential energy threshold of 
8.5 keV, are 2.8, 4.1, 4.7 and 4.9 nm. These calculations are repeated 
for the experiment presented here (Xe ions at 5 MeV) and 
0.004 MeV (slow HCI case^). The hereby obtained values of d/F 
are plotted in Fig. 5b against the experimentally measured Ep thresh- 
olds (5.5, 8.5 and 11.2 keV) with d as the parameter and are fitted 
linearly. Using the experimental value of d/F = 4.1 nm, the depths of 
interaction are 1.7, 2.0 and 2.3 nm, leading to F= 0.41, 0.49 and 0.56 
for Ep = 5.5, 8.5 and 11.2 keV respectively. It can be seen that the 
depth of interaction and the fraction of deposited potential energy 
increase with Ep within the i-TS model. Comparing this data to the 
values predicted by LemelP^, d is in the same range of depth but F is 
larger. 

To conclude, the required electronic excitation to induce a nano- 
hillock can be either deposited in the form of electronic energy loss by 
SHI or via potential energy by slow HCI or both, as in our case. We 
have shown that there is a clear quantitative linear relation between 
the increase of Sg and the decrease of the deposited Ep, leading to an 
additive action of the two mechanisms of energy deposition. This 
relation was quantified in a way that a link between the fraction of 
potential energy and the depth in which this energy is deposited can 
be estabUshed for CaF2: d (nm) = 4.1 X F. The inelastic thermal spike 
model applied to describe material modification by Sg supplied by 
SHI has been extended to the potential energy in three-dimensional 
geometry. Despite the number of parameters that are unknown, it is 
possible to predict values of the depth and the fraction of deposited 
potential energy. Measurements of the depth of the molten phase will 
be a crucial test of the inelastic thermal spike model applied to 
describe hillock formation by potential energy. 

Methods 

Preparation of the samples and AFM measurements. CaF2 (111) crystals (Korth 
Crystal Company) were freshly cleaved in air before placing them in the high vacuum 
irradiation chamber. After irradiation, samples were imaged under ambient 
conditions by Atomic Force Microscopy (AFM) in tapping mode using either a 
Nanoscopelll (DI) in Lanzhou or a Cypher AFM (Asylum Research) in Vienna. 
Image analysis was done with the Gwyddion software^^. Using statistical methods, 
this code allows determining the surface roughness and the number of hillocks per 
unit area for each sample. Freshly cleaved surfaces of CaF2 are flat with a mean 
roughness of 0.07 nm and no changes are observed when inspected after 5 and 16 
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Figure 5 | (a) Maximum energy deposited on the atoms (eV/atom) calculated with the i-TS model for Xe ions at 3 MeV versus the deposited 
potential energy (F*Ep) . The energy per atom at F*Ep = 0 corresponds to the electronic energy loss for Xe ions at 3 MeV. The deposited potential energy is 
given by a different cylinder volume with a radius of 1 nm for four different values of the depth (d = 1, 2, 3 and 4 nm). The dotted horizontal line 
is the energy per atom necessary to melt CaF2. (b) The calculated value of d/F when reaching the melting energy versus Ep for irradiations performed at 5, 3 
and 0.004 MeV is compared to the experimental ratio d/F = 4.1 nm (dotted line). 



days. Consequently, within the experimental errors, the irradiated surface keeps the 
same roughness over time. 

Medium energy irradiations. Irradiations were carried out on the 320 kV ECR 
platform for highly charged ions physics research at IMP ( Lanzhou) For each 
sample, the fluence was about 5 X 10^° Xe'^'^/cm^. CaF2 surfaces have been irradiated 
at two kinetic energies (3 and 5 MeV) by Xe'^"^ with charge states q between 19 and 26. 
All projectile charge states used were therefore below the charge state threshold for 
hillock formation reported to be Xe^*^"^ for very low kinetic energ)^ (0.004 MeV)^. In 
addition, irradiations with Xe^°"^ ions have been performed at three kinetic energies 
(0.54, 3 and 5 MeV). 

High energy irradiations. Swift heavy Xe ions were delivered by the IRRSUD beam 
line of the GANIL facility in Caen^*^. The initial beam energy was 92 MeV and using 
Al foil degraders, Xe ion energies ranging between 12 and 58 MeV could also be 
realized. The fluences applied were 5 X 10^ Xe/cm^ and 1.5 X 10^° Xe/cm^. 
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